Surface energy and stability of stress-driven discommensurate surface structures 
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A method is presented to obtain ab initio upper and lower bounds to surface energies of stress-driven 
discommensurate surface structures, possibly non-periodic or exhibiting very large unit cells. The 
instability of the stressed, commensurate parent of the discommensurate structure sets an upper 
bound to its surface energy; a lower bound is defined by the surface energy of an ideally commen- 
surate but laterally strained hypothetical surface system. The surface energies of the phases of the 
Si(lll):Ga and Ge(lll):Ga systems and the energies of the discommensurations are determined 
within ±0.2 eV. 
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In the present paper we show how surface energies of 
discommensurate surface reconstructions, which exhibit 
extraordinary large unit cells or are not periodic at all, 
can be determined with an accuracy of a few tenths of 
an eV using ab initio calculations of much simpler struc- 
tures. The method is applied to the Si(lll):Ga and 
Ge(lll):Ga systems, which display several discommen- 
surate reconstructions. 

Reconstructed surfaces can be classified according t< 
their driving force, which is most frequently chemical 
When the chemistry is satisfied, surfaces can still be un- 
stable towards rearrangements on a larger scale, driven by 
other factors like enthropyn or surface stress.ETU Stress- 
driven discommensurate surface-structures can be found 
for clean surfaces, like Au(lll),El and adsorbate systems. 
A rich phenomenology has been found .-for metals ad- 
sorbed on (111) surfaces of SiQ and Gea for coverages 
in the monolayer range. 

An instability towards discommensuration arises when 
the parent structure (the structure prior to discommen- 
suration) is-.under a sufficient compressive (tensile) sur- 
face stress.u A domain superstructure appears at the 
surface. Within domains the surface partially relaxes 
its stress by expanding (contracting) the surface lat- 
tice. The structure within domains resembles that of 
the parent structure whereas the domain walls, which 
allow the adjustment to the subsurface periodicity, rep- 
resent defects. They are also considered as excitations if 
the discommensurability occurs during a commensurate- 
incommensurate phase transitionfl 

The discommensurate structures of adsorbates on 
semiconductor surfaces studied here display several char- 
acteristic differences with respect to the better known 
physisorbed, discommensurate systems, like noble gases 
physisorbed on graphitejj which-are understood in terms 
of the Frenkel-Kontorova modeljl^l or with respect to the 
case of Au(lll), which is also understood in those, terms 
and considerations of continuum elasticity theoryfl'El The 
main differences are: (i) The adsorbate chemisorbs medi- 
ated by covalent, i.e., directional bonds, which represents 
a dramatic change of the energy scale, (ii) The chemistry 



changes drastically at the domain walls, (iii) The coex- 
istence of two discommensurate phases at certain cover- 
ages suggests a first order phase transition as function of 
coverageO [see Fig. 1 (a)]. These differences indicate the 
need for a new qualitative understanding. As a first step, 
we analyze the energetics of these complicated systems. 

First-principles studies of reconstructed surfaces are 
always computationally demanding. However, the sit- 
uation for discommensurate structures is worse since the 
surface periodicity is lost (or the unit cell is very large) 
and the geometry and bonding at domain walls may not 
be well characterized. Here we show how to obtain sur- 
face energies of stress-driven discommensurate surface 
structures, using only information obtained from com- 
mensurate structures of the system. We demonstrate 
this for the Si(lll):Ga and Ge(lll):Ga systems. Each 
system shows two stable discommensurate phases, 7 and 
(3, very close in saturation coverages but with quite dif- 
ferent structures [cf. Fig. 1 (a)]. Both systems are very 
similar, i-e.,m)<si is very similar to 7g C7 and /3gi very sim- 
ilar to /3gg -Qu However, they differ substantially in the 
fact that a commensurate -s/3 x v / 3i?30° phase (\/3 here- 
after) is stable for Si(lll):Ga but not for Ge(lll):Ga. 
This rich scenario makes these systems very interesting 
on their own, but also very illustrative of the uscfulncs 
of our approach. 

The first-principles methoyd-|Used here is based on the 
local density approximation,£3 and is described in detail 
elsewhere.Q A repeated slab geometry, 12 layers thick, six 
layers of vacuum, approximates the semi-infinite geome- 
try of the surface. fT-he surface energy (per surface unit 
a£3 



cell fi) is defined as 



E s = (E s iab — N su b s E i 



bulk 
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l )/2, 



(1) 



where E s i a i, is the energy of the slab system per slab 
unit cell, E^^ s is the energy per atom of the substrate 
(Si or Ge) in its bulk form, E^° m is the atomic energy 
of the adsorbate species (Ga), and N su b s and N a d are 
the number of atoms of substrate and adsorbate species, 
respectively, per slab unit cell. The factor 1/2 reflects 
the presence of two surfaces per slab. 
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A stress-driven discommensurate structure is stabilized 
(versus its parent) by the energy which is released as 
the surface is allowed to relax laterally within domains. 
However, the energy gain is reduced by the discommen- 
suration energy (Edi S ), which is always positive. Two 
terms contribute to E d i S : (a) The energy required for 
creating the domain walls, which represent defects from 
a chemical point of view, since they disturb the other- 
wise preferred atomic arrangement at the surface.E-3 (b) 
The energy which is necessary to allow for the local shears 
that appear for a strained surface on an unstrained bulk.cl 
While the exact values of these energy contributions are 
not accessible by todays ab initio methods, we can cal- 
culate upper and lower bounds of the surface energy of 
such discommensurate systems. We start with the latter. 

When the whole parent system, bulk included (the slab 
in this case), is strained laterally by e [Fig. 1 (b)], the 
energy of the system increases, but the surface energy 
EP ar (e) decreases towards a minimum [Fig. 1 (c)] at a 
strain value very close to the one observed experimentally 
within domains, e d i S - In the j3 phases half of the domains 
shcwrji stacking fault at the surface, the other half do 
not .LI B In this case EP ar (edis) is half of the sum of the 
surface energies of both strained parent structures (with 
and without stacking fault). Energies are obtained from 
first principles by relaxing the geometry up to the fifth 
layer for each value of the strain. Era Note that E^ l b k s in 
Eq. (0) is the energy of the substrate bulk under the 
same strain. Since the energies shown in Fig. 1 (c) refer 
to the surface unit cell of the strained system O' they 
have to be corrected by the ratio of areas of fl and f2', 
giving E l ° w = EP ar (e dls )/(l + e dls ) 2 . 

The lower bound is more precisely adjusted by tak- 
ing into account the small change of adsorbate cover- 
age at the domain walls 59 w (negative in the case of 
light walls like for Si:Ga and Ge:Ga). This is done by 
adding the term EP ar (edi S )S9 w /9 par , where 9 par is the 
Ga coverage of the parent structure. Since 1/(1 + e^js) 2 
= (@dis — 30w)/@par, where 9 dis is the Ga coverage of 
the discommensurate structure, we obtain the following 
lower bound: 

E l r = El ar {e dis )^. (2) 

"par 

The energy of a discommensurate phase \ 1S then ex- 
pressed as 

E x = E l r* + El s . (3) 

For Si: Ga! the coveragesO are # 7 = 0.8 monloayers 
(ML) and 9p = 0.3 ML, and 9 1 = 0.7 ML and O = 
0.8 ML for Ge:Ga.Q For the strain we use ^d^y^ 10%, 
similar to the experimentally observed valueurEl and in 
accordance with the minima in Fig. 1 (c). 

Upper bounds are obtained by an analysis of the sur- 
face energy as a function of adsorbate coverage, using the 
fact that the unstrained, commensurate parent structures 



are not observed experimentally. The stability of the-dif- 
ferent phases is shown by the convex-hull criterion£3 in 
the graph of energy versus coverage the lines joining the 
points representing stable, saturated surface phases de- 
fine a convex curve. This treatment is strictly equiva- 
lent to the analysis based on the grand-canonical frep. 
energy as a function of adsorbate chemical potentialjLZl 
both at zero temperature. Required for this analysis and 
observed experimentally are well defined saturation cov- 
erages for all the phases involved. 

The only variable-coverage Ga phase that grows on 
the (111) surfaces of Si and Ge is bulk Ga in the form of 
clusters. Assuming bulk-like, large clusters, their inter- 
nal energy per Ga atom, /iGa, can be approximated by 
the corresponding value in bulk Ga, /iot^.O This quan- 
tity has been calculated for the stable (orthorhombic) 
bulk Ga phase using the same technique and approxima- 
tions as for the surface calculations, giving = —3.45 
eV/atom (fJ,Q° m = 0). In the plot of surface energy ver- 
sus coverage, the cluster phase appears as a straight line 
with slope equal to /iQa fe ■ The line starts at the point as- 
sociated with the stable adsorbate phase of highest cov- 
erage, since the surface energy as defined in Eq. [l] is 
now the one corresponding to the adsorbate phase plus 
the cluster energy which is equal to HQg k (0tot — ad ). It 
should be noted that, depending on the growth condi- 
tions, Ga forms small clusters for coverages up to a few 
ML which may increase the effective /j,q" s (a kinetic ef- 
fect). 

To fullfil the convex-hull criterion, the curve which 
joins the points corresponding to saturation phases and 
ends in a line of slope iJq^ has to be convex. It is more 
convenient to use in the diagrams the reduced energy 

£ s = E s - M&f Ga . (4) 

The stability curve remains convex, but the cluster line 
has now zero slope (inset in Fig. 2). The-energiesJor the 
clean, reconstructed (111) surfaces of Silla and Geta, and 
\[2> phases foK-SiEj and for GeEJ have been taken from 
the literature. El 

In Fig. 2 convex-hull diagrams are shown for Si:Ga and 
Ge:Ga. The bars indicate the energy range of stability for 
the 7 and /3 phases, without taking into account the ex- 
perimentally determined (in)stability of the \/3 phases. 
The resulting surface energies are E] = —1.90 ± 0.32 
eV/fl and E? = -2.25 ± 0.32 eV/fl for Si:Ga, and 
E] = -1.71 ± 0.23 eV/fl and E% = -2.05 ± 0.22 eV/fl 
for Ge:Ga. The ± gives the margins of stability for these 
systems which are larger than typical errors of full first- 
principles calculations. However, the margins diminish 
with increasing domain size, since the strain is smaller. 
This reduces the difference in energy between parent and 
strained surfaces, the local shears, and the domain-wall 
energy per surface atom (lower density of domain walls) . 
Thus, in particular for very large domains, where full 
first-principles calculations become impractical, the ap- 
plicability and accuracy of our method improve. Fur- 
thermore, the relative energies of the discommensurate 
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phases are more precisely defined by the convex hull cri- 
terion: has to be within [£] — S,£]], and vice- versa, 
£] has to be within + S\, with S < 0.08 eV/il in 

this case. 

We obtain important information about the discom- 
mensuration energy Edi S which is the energy increase of 
the discommensurate phase compared to the lower bound 
(cf. Eq. g). As indicated in Fig. 2, Ej is > E dis 
by at least 0.10 eV/tt for Si and 0.14 for Ge [at most, 

E li S > E L h y (°- 10 + s ) eY / n for Si > and b y (°- 14 + S ) 

eV/fi for Ge]. This means energy differences of the or- 
der of 5 to 10 eV per domain. Taking into account that 
the contribution of the local shears to these energies is 
expected to be larger for j3 than for 7 (larger domains), 
these results indicate that we have chemically a better 
bonding situation at the domain walls for (3 than for 7, 
both for Si and Ge. This agrees with the experimetal 
findings of scanning tunneling microscopy measurements 
for Ge:Ga, showing saturated bonds, i.e. little density 
of states p. the band gap region at the walls for (3 but 
not for 7O This is actually the factor stabilizing the /? 
phase, a phase which is higher in energy within domains 
due to the presence of the stacking fault [cf. Fig. 1 (c)]. 

The problem of the (in)slability of the V3 phases is also 
addressed in this analysis. El Figure 2 shows that the -\/3 
phase in Si:Ga is stable, and that in Ge:Ga it is probably 
unstable (stable only if the energies for 7 and /3 would be 
at about their upper limits). Ell Experimentally, the a/3 
phase has been shown to be stable for Si but not for Ge, 
though it can be prepared for Ge at low temperature on 
cleaved (111) surf aces. c3 

The surface energies of the discommensurate phases 
can be restricted to a narrower range if we use the in- 
formation about the (in)stability of the y3 phases as 
inputEil (insets in Figure 2). The new surface energies 
are then E] = -2.05±0.18 eV/fl and = -2.39±0.18 
eV/n for Si:Ga, and E] = -1.76 ± 0.18 eV/tt and 
E% = -2.10 ± 0.18 eV/O for Ge:Ga with the ± giving 
again the margins. 

For Ge:Ga and approximately 1.5 ML Ga coverage, a 
third discommensurate structure, (3h 3 , has been charac^ 
terized experimentally and observed to be metastable.u 
It consists of Ga adatoms covering the (3 phase at H3 
positions. The lower bound is calculated as for the other 
structures and is found to be £° W '^ H 3 — q.84 eV/fi. It 
is just on the limit of stability (see Fig. 2). The surface 
strain present in the underlying (3 structure is an essential 
ingredient for the appearance of /3h 3 versus other possi- 
ble adlayer structures (see Fig. 8 in Ref. ||) . Its existence 
can be understood in terms of an energy barrier toward 
the formation of large Ga clusters (see above) and may 
furthermore suggest that the energies for 7 and (3 are 
closer to their upper than to their lower bounds. As dis- 
cussed in detail earlier,B we believe, that this structure 
mediates an easy transition from 7 to (3. 

In summary, we have shown how to assess the sur- 
face energy of discommensurate surface phases. Based on 



first-principles calculations of much simpler, commensu- 
rate structures and on stability considerations we obtain 
upper and lower bounds. We have demonstrated this ap- 
proach for Si:Ga and Ge:Ga. The (in)stability of the a/3 
phases has been also addressed and values for the ener- 
gies of the discommensurations have been obtained. This 
information is particularly important for understanding 
the phase diagram of these systems. 
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FIG. 1. (a) Scanning tunneling microscopy image of the 7 
(left) and 3 (right) phases coexisting on Ge(lll):Ga (taken 
from Ref. Illh . (b) Commensurate, strained commensurate, 
and discommensurate structures in side view, (c) Surface 
energy per strained surface unit cell (f2') versus lateral strain 
for Ge(lll):Ga-(lx 1) with Ga substituting the surface layer, 
with (full circles) and without (empty circles) stacking fault 
at the surface. D.W. means domain wall. 

FIG. 2. Stability diagrams for Si:Ga and Ge:Ga. Diamonds 
indicate unstrained commensurate structures, filled circles in- 
dicate lower bounds for 7 and B structures (and for 3h 3 in 
Ge:Ga), bars indicate their ranges of stability, and squares in- 
dicate \/3 structures. The insets show the diagrams in terms 
of the reduced energy (cf. Eq. 4). Arrows indicate ranges of 
stability if exploiting the knowledge about the (in)stability of 
the phases. 



4 



